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Abstract 

Solutions, in terms of special functions, of all wave equations Uxx — 
Utt = V{x)u{t,x), characterised by eight inequivalent time independent 
potentials and by variables separation, have been found. The real val- 
ueness and the properties of the solutions produced by computer algebra 
programs are not always manifest and in this work we provide ready to use 
solutions. We discuss especially the potential (mi + m2 sinh a;) cosh"'^ x. 
Such potential approximates the Schwarzschild black hole potential for 
even parity and its use for determining black holes quasi-normal modes is 
hinted to. 

PACS: 02.30Jr 04.70Bw 



1 Introduction 

In the frame of our research work on analytic solutions of black holes differential 
equations we have determined the analytic solutions of all wave equations 
with time independent potential V{x): 

mx (1) 

mx~^ (2) 

TO sin~^ X (3) 

TO sinh~^ X (4) 
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m cosh~^ X (5) 

m exp X (6) 

(mi + m2 sin x) cos~^ x (7) 

(toi + TO2 sinh x) cosh~^ x (8) 

(rrii + m2 cosh x) sinh~^ x (9) 

rriie^ + TO2e^'' (10) 

rrii + m2X~'^ (11) 

m (12) 

for the wave equation: 

Uxx - utt = V{x)u{t, x) (13) 

These potentials characterize the wave equation by variable separation . They 
may be reduced to eight irreducible forms. The potential lIHll is equivalent to 
potential l|12|l with the change of variables: 



exp (I) cosh ^ t' = exp sinh ^ 



while potentials (j3l4l5|l are equivalent to potential lO with the following change 
of variables respectively: 

= tan ^ + tan 77 i' = tan ^ — tan 77 (14) 

a;' = tanh ^ + tanh 77 i' = tanh ^ — tanh 77 (15) 

a;' = coth ^ + tanh 77 i' = coth ^ — tanh 77 (16) 



where 



The general form of the solution with separated variables of ea. p3|l is: 

u{t,x)=(t)l{LOl)(j)2{L02) (17) 

where wi — uJi{t,x), UJ2 — uj2{t,x), and 4>i{uji) and (^2(^2) are arbitrary solu- 
tions of the separated ordinary differential equation: 

J2 J, 

^ [{c^ + g^{uj^)](l), (18) 



duj 



Ci being the separation constant and i = t,x. 

The wave equations separate in several coordinate systems, among which 
the commonest and simplest poses lui = t, uj2 = x, gt — and = V{x). With 
the support of Maple 7 software, we have obtained the general solutions of all 
differential equations, analysed their properties and proved their real valueness. 
The solution for the time dependent function equation: 

^-.0,^0 (19) 
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is: 

= kie'^' + k2e-'^' (20) 

The space dependent solutions in the appendix have been obtained by simph- 
fying the results returned by the odesolve function of Maple 7. The solutions 
corresponding to potentials H7I8I9I) require an additional effort to obtain a read- 
able form^. It is the latter we directly show in the appendix. The applicability, 
the properties, and the real valueness of the solution corresponding to potential 
((HJ have been considered in detail in the next section. 



2 The Regge— Wheeler— Zerilli equation 

Regge and Wheeler T proved the stability of the Schwarzschild black hole in 
vacuum for axial perturbations, while Zerilli 5 found the equation for polar 
perturbations. The latter is written in terms of the wave function ^E*;, for each 
/-pole component : 

- - Vi{r)Mr,t) = (21) 

where 

r*=r + 2Afln(^-l) (22) 

is the tortoise coordinate and the potential Vi{r) is: 

, ( 2Af\ 2A2(A + l)r3 + 6A2Mr2 + 18AM2r + 18M3 
^("^^l^-— j r3(Ar + 3Af)^ ^'^^ 

while A = i(; - 1)(/ + 2). The Zerilh potential (ESJ may be approximated in 
a selected radial coordinate domain, including the maximum, by the potential 
Blome and Mashhoon^ have used the Eckart potential [7j: 

Y = Vbe^'" - yo{tanh[a(x - Xq) + /i] - tanh [lY cosh^ p. (24) 

while Ferrari and Mashhoon|Sj, Beyer 9^ have used the Poschl-Teller potentialjlO): 

V = ^ (25) 

cosh a{x ~ xq) 

for derivation of the QNM (quasi-normal modes) of a black hole. The ground 
state plus the first few excited states can be approximated by the bound states of 
the inverted potential. We note that ^ well reproduces the Zerilli potential and 
investigations on quasi-normal modes are imdergoing|ll|. The Zerilli potential 
has not allowed any analytic determination of the QNM^. The polar potential 
is thus substituted by: 

Vl[x) = A[mi + TO2 sinh(/cr*)] cosh"^(A;r*) (26) 



^Ref. 13 has been used for special function properties. 

■^Black holes perturbations equations do not admit exact solutions, apart of approximate 
solutions for portions of the frequency domain, e.g. 1121 — 118^ . or post— Newtonian expansions 
in the weak field and slow motion, e.g. |19| . 
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Figure 1: The Zerilli potential, the potential the Eckart and the Poschl- 
Teller potentials for I = 2 (quadrupole), starting from above. On the y axis the 
potential in (1/M^) units; on the x axis the radial coordinate in r/M units. 

where A,mi,m2,k are parameters depending on I and M, for proper curve 
fitting. 

From fig. 1, we evince that the original black hole Zerilli potential is best 
replaced by potential ^ for r/M < 10 and after by the Eckart potential. 
The general form of the solution is: 



The d'Alembert equation [df — d^* — V {r*)]'i/ 1 — separates into four coordinate 
systems, among which the following is the only with an explicit relation for tof. 

ujr' =r* LOt^t fr* = A[mi + TO2 sinh(fcr*)] cosh"^(fcr*) /t = (29) 

In order to write the solution corresponding to the potential JHJ in a more 
suitable way, we assume Cx <Q and we introduce new real parameters k, x, fj,i, 
fi2 and a such that: 

TOi = 1/4 — + /i2, 7712 = 2/ii/i2, A = k'^ , Cx = —a'^k'^, sinh(fcr*) = x 

Consequently: 



where ipHi i — r* ,t are solutions of the o.d.e.: 



(27) 




(28) 



as = 2(^1+1/^2) 



a ~ 1/2 + fii ~ i(7 

b = 1/2 + fii+ia 

c = 1 + fli + ifl2 

z = (X-ix)/2 
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Accordingly, the potential (|23|l can be rewritten as follows: 

where the function w{z) has to be a solution of the hypergeometric equation: 

z{z - 1) + [z(l + 6 + a) - c)]— ^ + a6u.(z) = (31) 

az^ clz 

and the general solution of which is: 

w{z) = kiF{a, b; c; z) + fc2Z^"'=F(a - c+ l,5-c+l;2-c;z) (32) 
Ea. H28|l has, therefore, the following expression: 

= [C0sh(fcr*)](^^+^/^) g-p.arctan[l/si„h(/c.*)]^ J^l-^Smhfcr*^ ^^^^ 

The solutions (|32|l must be real valued functions. This is possible since: 

a = h, b = a, c — a + b+1 — c, z — 1 — z (34) 

so that: 



F{a, b; c; z) — F(a, b; c; z) ~ Fib, a; a + + 1 — c; 1 — z) (35) 

and, moreover, for one of the identities between hypergeometric functions we 
know that: 

F(6, a; a + + 1 — c; 1 — z) = <JiF{a,b]c; z) (36) 

+(72Z^^''F{a - c + 1, - c + 1; 2 - c; z) 

where we set: 

^ r(c)r(i - c) ^ r(c)r(i - c) 

' r(a-c + i)r(6-c + i)' r(a)r(6) ^ ^ 

Finally, for the potential Q the space dependent solution is: 

0, = 2(c/2-l/4)(^_^)[(a+&-c)/2+l/4] 

fciF(a, 6; c; z) + fc2z(i"'=)F(l + a - c, 1 + 6 - c; 2 - c; z)l (38) 



where 



a = (ag + a4)/4 + /3 + 1/2 

6 = (ag + a4)/4 - /3 + 1/2 

c = 1 + 03/2 

z = (1 ± i sinha;)/2 



with ag = ±\/l — 4mi ± 4ir7i2, a4 = ±\/l ^ 4?7ii =p 4iTO2 and /? = ±y^. 



3 Conclusions 



Solutions, in terms of special functions, of all wave equations u^x — uu = 
V{x)u, characterized by eight inequivalent time independent potentials and 
by variable separation, have been found. One among the potentials, (mi + 
7Ti2 sinha;) cosh~^ X, has a shape similar to Schwarzschild black hole potential 
for polar perturbations. 
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5 Appendix 

In this section we write the space dependent solutions for all potentials, obtained 
after some computation. For the potential the solution can be written in 
terms of the Airy functions: 



^ ^ fc^Ai ( ^^4^ ] + fc2Bi ' ''''' 



or in terms of Bessel K and I functions: 



= VCx + mx <j fciK ( -, - 



,2/3 



3 \m\ 
1 2{cx + mxf/'^'' 
3' 3|mj 



1 2(c^ + mxf/'^ 



3' 



3|m| 



(39) 



(40) 



For the potential Q: 



KlJ I .X\/-C2 



(41) 



where J and Y are Bessel functions. 
For the potential ©: 



V: 



sm X 



1 1 



fciP ( iy/c^- ^, ^Vl + 4m, COS a; 



1 1 



-hQ ( i^/c^- -, -Vl + 4m,cosa; 



(42) 



where P{v,fi,z) = P'^{z) and Q,{v,^,z) — ^^{z) are associated Legendre func- 
tions of the first and second kind, respectively. 
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For the potential 



V sinh : 



1 1 



kiP ( a/cT- -, -Vl + 4m, cosh a; 



1 1 



-fc2Q ( \/c^- 2' + 4m, cosh a; 



(43) 



For the potential 

(hx ~ V cosh X 



1 1 



fciP ( -,-vT^^4m, z|sinhx| 



1 1 



+fc2Q ( \/c^- -, -Vl - 4m,i |sinha;| 



(44) 



For the potential ©: 

= A:i J (273;, e"/2^ + fcsY (2^^;, 2^^ e"/^^ (45) 

For the potential {Tj): 

Cj)x = ^(c/2^1/4)(^_^)[(a+b-c)/2+l/4] 

fciF(a, 6; c; z) + fc2z(^"'^)F(l + a - c, 1 + 6 - c; 2 - c; z) (46) 

where F indicates Gauss hypergeometric function 2F1 and the parameters a, h, 
c and z are defined by: 

a = (ai + a2)/4 + /3 + 1/2 
h = (ai + a2)/4 - /3 + 1/2 
c = 1 + 0^2/2 
z = (l + sina;)/2 



with ai = ±-\/l + 4mi + 4m2, a2 = ±-\/l + 4mi — 4m2 and /3 = ±\/"^c^ 
For the potential ©: 



(z-1) 



[(a+&-c)/2+l/4] 



fciF(a, 5; c; z) + fc2z(^"'^)F(l + 6 - c, 1 + a - c; 2 - c; z) (47) 



where 



a = (as + q:6)/4 + /3 + 1/2 

& = (as + a6)/4 - /3 + 1/2 

c = 1 + a6/2 

z — (1 + cosha;)/2 



with as = ±^/l + toi + m2, ag = ±-\/l + m-i ^ m-2 and /? — ±^/cZ. 
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For the potential ((Tn|l : 

-a;/2 



e 



fciM ( --^^7^, -2V^ ) + fcaW (-77^, V^, -2^^ 

(48) 

where W and M are Whittaker's functions. 

For the potential the space dependent solution is equal to the solution 
for the potential jSJ, while for the potential (|12|1 . the solution is trivial: 

(j)^ = i^ie^^ + ifae-v^^ (49) 
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